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Abstract. The paper deals with several properties of the Oort cloud of comets. Sun, Galaxy (and 
("V^ ' Jupiter) gravitationally act on the comets. New physical model of galactic tide is considered. 

Cn , The main results can be summarized as follows: 

^.,^' 1. Mass of the Oort cloud of comets is less than 1 mass of the Earth (Me), probably not greater 

C^ I than 1/2 Me- 

2. Theoretical number of long-period comets with perihelion distance q < 5 AU is about 50- 
times greater than the conventional approach yields. Gravity of Jupiter was taken into account 

^\ , in finding this result. 

3. Semi-major axis a and period of oscillations P of eccentricity (and other orbital elements) 
are related as a'^ P = 1 in natural units for a moving Solar System in the Galaxy. The natural 
unit for time is the orbital period of the Solar System revolution around the galactic center and 
the natural unit for measuring the semi-major axis is its maximum value for the half-radius of 
the Solar System corresponding to the half-radius of the Oort cloud. The relation holds for the 

rS^ . cases when comets approach the inner part of the Solar System, e.g., perihelion distances are 

Ph' less than « 100 AU. 

4. The minimum value of semi-major axis for the Oort cloud is amin -^ 1 x 10* AU. This 
^ ■ condition was obtained both from the numerical results on cometary evolution under the action 

"v^ ' of the galactic tides and from the observational distribution of long-period comets. If the density 

d , function of semi-major axis is approximated by proportionality a", then a is — 1/2, approxi- 

mately. 

5. The magnitude of the change in perihelion distance per orbit, Aq, of a comet due to galactic 
►-^ , tides is a strong function of semi-major axis a, proportional to a^'^^. 
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1. Introduction 

Intense modeling of the Oort cloud of comets started practically immediately after the Oort's paper 
(Oort 1950). As for the curent status of the ideas about the Oort cloud we refer to review papers, e. 
g., Dones et al. (2004), Levison and Jones (2007). 
. ^_, ^ New physical access to the theoretical modeling of the Oort cloud of comets was suggested by 

^^ ■ Klacka (2009a, 2009b). It is based on a new approach to treating the effect of galactic tides, including 

H I improvement of the model of the Galaxy. Its significance was discussed by Komar et al. (2009). 

Analytical approach to secular evolution of orbital elements was discussed by Pastor et al. (2009). 
The aim of this paper is to present some characteristics obtained from the new model. The new 
model is represented by Eqs. (26)-(27) in Klacka (2009a) and denoted as the Model II in Komar et al. 
(2009). The last paper showed the importance of the new model in comparison with the conventional 
models. Thus, we will not make any calculation based on the conventional models. Instead of that 
we will concentrate on calculations based on the new modeling of galactic tides. The results of the 
calculations will be compared with the results obtained by other authors isung conventional models 
of galactic tides. 

Sec. 2 deals with distribution of cometary inclinations and it sheds some light also on the value of 
the exponent a describing distribution of semi-major axes of comets in the Oort cloud. Sec. 3 shows 
that the magnitude of the change in perihelion distance per orbit of a comet due to galactic tides 
is a strong function of semi-major axis a, proportional to a^'^^. Sec. 4 relates the semi-major axis of 
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a comet and period of oscillation in other orbital elements due to galactic tides. The found relation 
reminds the third Kepler's law, also in the case when natural units of length and time are used. Sees. 
5 and 6 discuss some theoretical access to cometary distributions in eccentricity, perihelion distance 
and semi-major axis. However, Sec. 6.3 presents also results obtained by numerical integration of 
equation of motion for comets under the action of the Sun and the Galaxy. Sec. 7 determines minimal 
perihelion distances Qmin and inclinations for the instant when qmin occurs. The section uses analytical 
approach to secular time derivatives of orbital elements when the initial cometary inclinations differ 
from 90°. The two sources of gravity, the Sun and the Galaxy, are used also in Sec. 8 dealing with 
the distribution of comets in inclination to ecliptic. Sec. 9 improves the most relevant previous results 
also for the case when the planet Jupiter gravitationally influences the motion of the comets. Finally, 
Sec. 10 deals with the distribution function in semi-major axis and the number of comets in the Oort 
cloud, including estimate of the mass of the Oort cloud. 

2. Distribution of cometary inclinations 

We are interested in distribution function of inclination of comets. The inclination is the angle between 
a reference plane and the instantaneous orbital plane of a comet. 

Probability that a comet has an inclination lying in the interval {i — di/2, i + di/2) is 

dp ~ h{i) dfl , 

dfl = 2 TT sini di , (1) 

where dfl is the solid angle and h(i) is a function, i e (0, tt). It is assumed that the normalization 
condition 



h{i) dn = 1 



4 7r 



2 TT / h{i) sin i di ^ 1 (2) 

Jo 

holds. The distribution function is 

H{i) = 2 TT / h{i') sini' di' , 

Jo 
H{n) = 1 . (3) 

The density function is 

, , dH 1 / \ 

huM = —r- = 2 vr hU) smz . 
di 

Hin) = 1 . (4) 

The mean value is given by relation 

(i) = 2 TT / i' h{i') sini' di' , (5) 

Jo 

The special case of the isotropic distribution is defined by h{i) = constant. Eqs. ©-dS]) reduce to 
H{i) = 2 (1 ^ '^°'^*) ' 
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hnii) = 
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hiz) E 


1 

' 47r ' 



(6) 
Equal numbers of prograde (i < 7r/2) and retrograde (i > 7r/2) orbits exist. 
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2.1. Results of orbital evolution - inclination to galactic equator 

We are dealing with distribution of comets in the Oort cloud in this subsection. We are interested in 
the distribution of inclination measured with respect to the galactic equatorial plane. 

We made detailed numerical calculations of the orbital evolution of comets under the action of 
gravity of the Sun and Galaxy (Klacka 2009a, 2009b, Komar et al. 2009). The results showed that 
initial inclination of a comet decreased to some value if Zi„ was smaller than 90° and increased to an 
another value if ii„ was greater than 90° . In reality the inclination changes on a time scale of millions 
of years. Fig. 1 depicts such kind of evolution of inclination of a comet moving initially on almost 
circular orbit in a distance of 5.0 x 10^ AU. We are interested in time average of the inclination in 
order to find results comparable with statements presented in literature (e.g., Duncan et al. 1987, 
Bailey 1983, Fernandez and Ip 1987, Fernandez 1992, Fernandez and Gallardo 1999). 
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Fig. 1. Time evolution of inclination of a comet moving initially on almost circular orbit in a distance 
of 5.0 X 10^ AU. Three intial values of inclination are considered. 

Table 1 presents the averaged time values of a comet with a given initial inclination. On the basis 
of the data presented in Table 1 we can conclude that the function h(i) is not a constant. E.g, the 
interval ii„ S (70°, 90°) produces the time averaged values i £ (60°, 90°). An approximation yields 
that for a mean value i[deg] w 60 + (90-60)/2 = 75 the following value holds: h (75°) w (90 - 70)/(90 
— 60)/ (4 tt) = 2 / 3 / (4 vr). Similarly, the interval ii„ € (25°, 70°) produces the time averaged values 
i e (25°, 60°). An approximation yields that for a mean value j[deg] sa 25 + (60— 25)/2 = 42.5 the 
following value holds: h (42.5°) « (70 - 25)/(60 - 25)/ (4 tt) w 4/3 / (4 n). Thus, we can conclude 
that /i(7r/4) « (4 / 3) (4 tt )^^. One can see the difference between the value and the result presented 
in Eqs. ©. The result /i(7r/4) « (4 / 3) (4 tt )-i holds for semi-major axis a = 5.0 x 10'' AU. The 
smaller a, the closer the result to that given in Eq. ©. 
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0.00 
10.00 
30.00 
50.00 
70.00 


0.00 

9.89 

29.31 

45.12 

60.28 


90.00 
110.00 
130.00 
150.00 
180.00 


89.47 
119.72 
134.88 
150.69 
180.00 



Table 1. Values of time averaged inclination i for a comet with a given initial inclination 
initially almost circular orbit and semi-major axis a — 5.0 x 10"* AU due to the galactic tide. 



for 



Detailed numerical calculations show that the function h defined by Eqs. ([T]) depends not only on 
the inclination i, but also on the semi-major axis a. The most relevant results are presented in Table 
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2. The function for a — 2.5 x 10* AU is depicted in Fig. 2, where also isotropic function is shown, for 
comparison. 
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{ 0, ^/2 ) 
0.71 X 7r/4 
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1.0 
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Table 2. Maxima of the function /i(i, a), hmax = h{imax, a), for various values of cometary semi-major 
axes a. The maximum of h holds for the angle imax- The interval i € ( 0, 7r/2 ) is considered. 




Fig. 2. Function of inclination with respect to the galactic equatorial plane. The solid line holds for 
the gravity of the Galaxy and the Sun and for the semi-major axis a = 2.5 x 10'' AU. The dotted line 
corresponds to the two-body problem. 



The function can be approximated by the following conditions: 
, 



dh(i,a) 
lim 



i-fO di 

/.(|,a)=0 

lim ^M!l^ = 

i^TT/2 di 



lim ^M^ = 



h{imax,a) = ma.xlh{i,a),i G (0, -)| , 



di 



(7) 



Moreover, 



h\^-^ + i,aj = hi^- - i,a 
^ G (0,|>. 



(8) 



On the basis of the values presented in Table 2 and Eqs. ([7]) we can make the following approxi- 
mation 



h{i, a) — ha{i) x {1 + i + Ci2 a i + C22 a i + C13 a i + C23 a i ) 



1, 



hail) = (2- 2J (-) « + (7^-3) [-) I- - ^ 



C12 



4.89184 



Klacka et al.: Properties of the Oort cloud 



C22 = - 0.80376 , 
ci3 = - 3.23854 , 
C23 ^ + 0.67758 . 

The function independent of semi-major axis is 
J a"h{i,a)da 



(9) 



/ a" da 

Q-min 

a e (-4,-2); - 3/2 



(10) 



We inserted also the values of the exponent a taken from Fernandez and Gallardo (1999), Duncan et 
al. (1987). We finally obtain 
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The quantity K = Kia,amin,amax) is determined by the normalization condition 



2 TT 



hii,a,armn,amax) Slli i di 



(11) 



(12) 



Having the function defined by Eqs. (jlip - p2p . we can find the distribution function defined by Eq. 
([3]) and the observational mean value (i) . It follows 



^ ' 2 



(13) 



Our detailed numerical calculations considering tidal effect of Galaxy show that amplitude of semi- 
major axis a rapidly increases with the value of a when a > 5 x 10'' AU. If a = 7.5 x 10"* AU, then 
the amplitude of oscillation is 5 x 10^ AU. Moreover, if a is comparable with 1 x 10^ AU, then an 
increase of semi-major axis exists. More correctly, there does not exist a constant value of semi-major 
axis for secular evolution, an amplitude is about 25 x 10"^ AU. Thus, the effect of galactic tide causes 
that stability of the Oort cloud exists only for a < Umax and Umax ~ 0.8 x 10^ AU. This seems to be 
consistent with observations, according to which Umax ~ 1 x 10^ AU and amin ~ 1 x 10^ AU (see 
also Fig. 1 in Fernandez 1992), although the value amin = 1-5 x 10** AU is also presented (Fernandez 
and Ip 1987). 

Fig. 3 depicts the function(s) ft,(z). Although the presented curves do not exhibit any significant 
difference, the real results are relevant. At first, if we take into account amin ^ 0, more correctly amin 
<C 1 x lO'' AU, then various values of a can be used and the function practically does not depend 
on the real value of a. The value amin — >■ is consistent with the data obtained from observational 
data on long-period comets (see Sec. 10), but it is not consistent with the conventional statements 
discussed in the previous paragraph. However, if we would like to use the values not fulfilling amin — > 
0, e.g., amin = 1-0 X 10'* AU or amin — 2.0 x lO"* AU, then the only acceptable value of a is — 1: the 
value a ^ — 1 would produce nonpositive function h. The value a = — 1 is not consistent with the 
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conventional statements on the Oort cloud that a — — 3/2 or a G ( — 4, — 2 ) (Duncan et al. 1987, 
Bailey 1983, Fernandez and Ip 1987, Fernandez 1992, Fernandez and Gallardo 1999). Moreover, the 
observational data suggest that a is closer to zero: a = 0.13 (see Eqs.[74]in Sec. 10), or a = — 0.55 
(see Eqs.ESlin Sec. 10). 





Fig. 3. The function of inclination with respect to the galactic equatorial plane h(i). Averaging over 
distribution in senii-niajor axis is done. The upper/left part of the figure is characterized by the values 
a,min — > 0, amax = 10 X 10'* AU (practically independent of a) . The lower/right part of the figure 
corresponds to the case amin = 1 x 10^ AU, amax = 10 x 10** AU and a = — 1.0. 



3. Timescales on which a comet's perihelion changes 

A formula for a timescales on which a comet's perihelion changes is presented by Levison and Dones 
(2007, p. 583). However, as it is discussed by Komar et al. (2009), the formula does not correspond 
to reality. It is important to find mathematical relation(s) which are based on physical approach. 

In order to obtain the cases when comets approach the inner part of the Solar System, e.g., 
perihelion distances are less than « 100 AU, we consider initial inclinations with respect to the 
galactic equatorial plane about 90 degrees. 



3.1. First approach 

Let us consider the change of perihelion distance per revolution of a comet around the Sun. The period 
of revolution of the comet is T [years] , the change of perihelion distance during the period T is Ag 
[AU] and semi-major axis of the comet is a [AU]. Fig. 4 holds for the cases when q obtains its minimum 
values during the time evolution q{t) for 4.5 x 10^ years. The solid line is an analytical approximation 
to numerical results obtained from evolution for the model of Klacka (2009a, Eqs. 26-27; Model II in 
Komar et al. 2009). We can conclude that 



logio ( ^ j 

A2 



Al + A2 logiQC 

- 35.95 ±0.89 , 
-F 6.75 ±0.20 . 



(14) 



3.2. Second approach 

We want to find a formula which is more similar to the mathematical formula presented by Levison 
and Dones (2007, p. 583). As it was already stressed, our formula must respect physical reality. We 
can use detailed numerical calculation for orbital evolution for our new model (Klacka 2009a, 2009b; 
Model II in Komar et al. 2009), or we can try to use the analytical approach to secular evolution of 
orbital elements (discussed by Pastor et al. 2009). The analytical approach is limited by the condition 
that semi-major axis of a comet must be less than about 1.5 x 10'' AU if the condition T x ( 1 / T^ + 
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Fig. 4. Change of perihelion distance Aq [AU] per revolution of a comet around the Sun. The period of 
revolution of the comet is T [years] and semi- major axis of the comet is a [AU] . Solid line corresponds 
to linear fit of the numerical solution (dashed curve) for the tidal effect of the Galaxy. 



1 / To) < 0.05, where T^ is the period of oscillations of the Sun with respect to the galactic equatorial 
plane and Tq is the period of revolution of the Sun with respect to the center of the Galaxy. 
Detailed solution of the equation of motion corresponding to Model II yields 



logio ^ = ^1 + ^2 logio a + As logio Q > 



Al 

A2 

As 



33.67 ±1.48 , 
6.29 ±0.31 , 
0.12 ±0.07 . 



(15) 



Again, as in the previous subsection, the period of revolution of the comet is T [years], the change 
of perihelion distance during the period T is Aq [AU] and semi-major axis and perihelion distance of 
the comet are a [AU] and q [AU]. Eqs. ([T5|) contain the same quantities as the relation presented by 
Levison and Dones (2007, p. 583). It is immediately seen that the values of A2 and ^3 significantly 
difi^er from the values presented by Levison and Dones. 



3.3. Comparison of the approaches 

We have found two relations, one of them is represented by Eqs. ([Ti)l . the other one by Eqs. p^ . None 
of them is consistent with the relation presented by Levison and Dones (2007, p. 583). Moreover, the 
relation represented by Eqs. (jlSp is characterized by much greater errors of the coefficients than it is 
in the case of Eqs. (fH)) . 

Secular approach to the orbital evolution treated by Pastor et al. (2009) also confirms that the 
relation presented by Levison and Dones (2007, p. 583) is not consistent with physics of the galactic 
tide. 

According to Eqs. (fH|) - P^ we see that the magnitude of the change in perihelion distance per 
orbit, Aq, of a comet due to galactic tides is a strong function of semi-major axis a, proportional 
to a*-^^, or considering some error, a^ , S e (7.5, 8.5). This is a new result, not consistent with the 
conventional value 5 = 3.5 (see, e.g., Dones et al. 2004, p. 155). 



4. Relation between semi-major axis and oscillation period 

Galactic tide causes oscillations of eccentricity (perihelion and aphelion distances, angular orbital 
elements) of a comet in the Oort cloud (see Fig. 7 in Komar et al. 2009). We are interested in a 
relation between the semi-major axis a and the oscillation period P of the comet. 

We are interested in the cases when comets approach the inner part of the Solar System, e.g., 
perihelion distances are less than w 100 AU. Thus, we consider the values of about 90 degrees for 
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the initial inclinations with respect to the galactic equatorial plane. Various values of other orbital 
elements are taken into account. 

Detailed numerical calculations for 1 x 10* AU < a < 8 x lO"* AU yield the result presented in 
Fig. 5, which yields 

logio a = A logio P + B , 

A = - 0.362 ±0.029 , 

B = + 7.846 ±0.254 , 

[a] = AU , [P] = yr . (16) 




log P [years] 



Fig. 5. Semi-major axis a as a function of oscillation period P of eccentricity of a comet from the 
Oort cloud. The oscillations are caused by the effect of galactic tide. Solid line corresponds to linear 
fit of the data obtained by numerical solution for the tidal effect of the Galaxy. The data for initial 
inclination 90 degrees are used (other values of inclinations may yield a little different values) . 



Since the period of oscillations slightly changes on the scales of billions of years, we can rewrite 
Eq. P^ into a more simple relation a^ P = constant. The value of the constant is determined by the 
least square method. We obtain 



a P — constant , 
constant = (6.93 ± 0.23) x lO^^AU^ yr , 
[a] = AU , [P] = yr . 

In natural units for the moving Solar System in the Galaxy we obtain 

a^ P = 1 . 



(17) 



(18) 



The natural unit for time is the orbital period of Solar System revolution around the galactic center, 
27r/{A — B) = 2.3 x 10* yrs. The natural unit for the semi-major axis is its maximum value for the 
half-radius of the Solar System which equals to the half-radius of the Oort cloud, 0.7 x 10^ AU. The 
value yields limiting cometary aphelion distances 1.4 x 10^ AU, half of the distance between the Sun 
and its nearest star. 

As a first approximation, a comet of the Oort cloud revolves around the Sun on an ellipse ful- 
filling the third Kepler's law c? jT'^ — constant{two body problem), where T is period of revolution 
of the comet. Moreover, eccentricity of the comet oscillates due to the effect of galactic tide and 
the period of oscillations of the eccentricity relates the semi-major axis of the comet as a'^ P = 
constant{galactic tide). 
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5. Probability that eccentricity is less than a given value 

We are interested in the probability P(e < cl) that eccentricity e of comets is less than a given value 

Orbital evolution of eccentricity due to the effect of galactic tides is a periodic function of time, 
approximately. We will use the following approximation: 



eit) 
e{t = 0) 



i |l - sin 
_ 1 

e^n - 2 

2 . , 

— arcsin (26^. 





\ "1 


5™1 


sin ( 2 '5»«) 


5 



1) 



Sin G (- i,+i) , e™ e (0,+l> , 



(19) 



for the time evolution of eccentricity e{t) with the period of oscillations P, as it is discussed in the 
previous section. 

5.1. Density and distribution functions of eccentricity 

Let an initial distribution of eccentricities of comets in the Oort cloud is represented by a density 
function /(ei„). We are interested in the density function /(e) at a given time t. We have e = e(t), 
e(t = 0) = Cin. Gravitation of the Sun and Galaxy will be considered. 
If F is the distribution function, then 



dF = f[e{t)] de , 

dF = fiCin) dCin . 

Knowing the inverse function e^ = g 

Let us rewrite Eq. p^ into the form 
1 



"^(e) to the function g{ein), Eqs. (1^ yield 



Ke^n) = 



l_2smf— j 



y/ein{l-ein) + cos f — j (2ei„ - 1) 



(20) 



(21) 



(22) 



Since this function is not monotonous, it's inverse function does not exist. However, it is possible 
to split the function g{ein) into two monotonous functions. There exist inverse functions for the two 
monotonous functions: 

1 + (2e - 1) cos ( -— 



and 



(e) = 


- v/e(l - e) 


. f27rt\ 


1 
+ 2 


em e 


(0,sm^(f)), 




[ 




. f27Tt\ 

sm 

\P J 


1 

+ 2 


(e) = 


+ Ve(l-e) 



(23) 



1 + (2e - 1) cos 



27ri 
~P" 



e sm^ 



,1) 



(24) 



Eqs. (EID, (ESD-dMl) yield 
1 



/(e) = f[g-'{e)] 



|cosx 



Xel\ 
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(2e— l)cosx sinx — 2-\/e(l — e) |sinx| sinx 
^l + [l + 8e(e-l)]sin2x- (2e- 1)2 + 4(2e - 1)^6(1 - e) |sinx| cosx 

Sin = .g"^(e) = - \/e(l-e) |sinx| + - [1 + (2e - l)cosx] , 

2 TT 

e,„ e (Csin^f^V , (25) 



and, 



f{e)^ f[g-'{e)] ^ 



cosx + X, 



e2\ 



(2e— l)cosx sinx + 2-y/e(l — e) sin x 

A„2 = 



l + [l + 8e(e-l)]sin2x- (2e- 1)2 - 4(2e - l)v/e(l - e) |sinx| cosx 



em = 5 ^(e) = + Ve(l-e)|sinx| + - [1 + (2e - 1) cosx] , 
2 -K 

e^n e (sin^ f^j , 1) . (26) 

The distribution function is 



F(e) = / /(e') de' , 



F{1) = 1 . (27) 

We may introduce that instead of Eqs. (QUI) an alternative approach can be used. Knowing the 
density function /(ei„) and e = g{ein), we may use a convolution for finding the density function /(e). 
The convolution yields 

/(e) = J /(e„) 5 [e - .g(e„0] de,„ . (28) 

5.1.1. Uniform distribution 

Let us consider, as an example, a uniform distribution of eccentricities. The density function is 

/(e„0 = 1 , e^e (0,1) . (29) 

The distribution function is F(ei„) = /q'" /(e^„) (ie^„ — e™, F{1) = 1. 
Using Eqs. ([25]) -(US]), we obtain 

/(e) ^ 



cosx + Xel\ 



(2e— l)cosx sinx ~ 2^e(l — e) |sinx| sinx 

Xei — = , 

Y/l + [l + 8e(e-l)]sin2x- (2e- 1)2 + 4(2e - 1) Ve(l - e) |sinx| cosx 

2 TT 

e„ G (0,sin2('^V , (30) 
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and, 

f{e) = f[g~He)] ^ 



cosx + X, 



e2\ 



(2e— l)cosx sinx + 2-\/e(l — e) sin x 
^e2 — I = , 

^Jl + [l + 8e{e~l)]sm^X~{'2e~l)^ - 4(2e - l)\/e(l - e) |sinx| cosx 

2 TT 

X = ^t, 
e^n e (sin^ f^j , 1) . (31) 

As an illustration we take t = P/2. Eqs. ([5IH) - ([5T|) yield /(e) = 1, and, the distribution function is 
F{e) = e: 

/(e) = 1, i = 0, 

F(e) = e , i = , 

/(e) = 1 , i = P/2 , 

F(e) = e , t = P/2 . (32) 

5.1.2. Distribution function F{ein) = ef^ 

We will consider, as an another example, F(ei„) = ef„. Motivation for this distribution function 
comes from Fernandez and Gallardo (1999), Hills (1981) and Jeans (1919). The corresponding density 
function is 

f{ein) = 2 dn , e„i e (0, 1) . (33) 

Using Eqs. ([25]) -([261), we obtain 



1+ 2e-l cosx - 2^e 1-e smx 
|cosx + Xei\ 
(2e-l)cosx sinx - 2^e(l - e) |sinx| sinx 
Y/l + [l + 8e(e-l)]sin2x- (2e- 1)2 + 4(2e - 1) Ve(l - e) |sinx| cosx 

2 TT 

X = —t, 



G (0,sin2(^)), (34) 



and, 



^ l + (2e-l)cosx + 2^e(l-e)|sinx| 
^^^' Icosx + Xe2| 

(2e — l)cosx sinx + 2^e(l — e) sin^ x 

^e2 = 



l + [l + 5e(e-l)] |sinx| sinx- (2e- 1)2 - 2(2e - l)v/e(l - e) |sinx| cosx 

2 TT 

e„, G <^"^'(^) '1) • (3^) 

As an illustration we take t — P/2. Eqs. (|34p - ([55|) yield /(e) = 2 ( 1 — e ), and, the distribution 
function is F{e) = e { 2 ~ e ): 

/(e) = e, i = 0. 



12 
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Fie) ^ - e^ 


, i-0, 


/(e) - 2 (1 


- e), t = P/2, 


F(e) - e (2 


- e) , t^P/2. 



The situation represented by Eqs. (1551) is illustrated in Fig. 6. 



(36) 
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Fig. 6. Density /(e) and distribution F{e) functions of eccentricity for two instants: i = and t = 
P/2, where P is the oscillation period of eccentricities (perihelion and aphelion distances, inclinations, 
...). The oscillations are caused by the effect of galactic tide. 



5.2. P{e < sl) for a random time 

Now, we want to take into account some other important facts. At first, we do not know at which time 
a given distribution /(ei„) holds. Then, we know from the previous section that there is a relation 
between the semi-major axis a and period of oscillations P. Thus, comets with various values of a are 
characterized with various values of P. Moreover, we do not know the exact value of the parameter 
a defining the distribution in semi-major axis, as it is presented in Sec. 2. All these facts lead us to 
the conclusion that we have to find an another possibility of obtaining the probability P(e < e^,) that 
eccentricity of comets e is less than a given value cl ■ 
Let us consider Eqs. (fTOj) . We can write 



P{e < cl) 

AtL 



P 

\tL 1 



tL 



ei 



^' 



Klacka et al.: Properties of the Oort cloud 



2 4 - - 5„ 



The last equation for cl determines the two values II i and t^ 2- Thus, we finally obtain 



Pie < eA 



1 



— arcsin (1 — 2 e^) 



13 

(37) 

(38) 



6. Probability that perihelion distance is less than a given value 

We are interested in the probability P{q < ql) that perihelion distance q of comets is less than a 
given value ql- We will assume that all comets have the same semi-major axis a. Sec. 6.4 will consider 
joint/cumulative distribution function consisting of two marginal distributions, one in a and another 
one in q. 

Let the orbital evolution of eccentricity due to the effect of galactic tides is given by Eqs. P^ . 
Since g = a ( 1 - e ), Eqs. (HH) yield 



1 



r. ■ f2TTt\ I (2-Kt 

+ 2sm I — 1 Vem(l - ei„) - cos I — 



(2e,„ - 1) 



(39) 



6.1. Distribution and density functions of perihelion distance 

We are interested in the probability P(q < q^) that a perihelion distance q is less than a given value 
qL- Since q — a{l — e) and semi- major axis a is a constant in secular evolution, we immediately 
obtain P{e < cl) = P{q > qL)- Thus, 

P{q < <zl) = 1 - P{e < bl) , 



F,{q) 



1 



F{e) 



P{q < qL) - F.iqL) , 

a 

q G (0,a) , 

where Fg is the distribution function of perihelion distance. 

The density function fq{q) of perihelion distance q is, on the basis of Eqs. (|40)) . 

fq{q) dq ^ - /(e) de , 

fM) = - /(e) , 
a 

e-l-'-, 
a 

q G (0,a) . 
Eqs. (HSI-dMl) and Eqs. (gT]) yield for the density function fq{q): 



(40) 



(41) 



fqiq) = - /(e) , 9G (0,a) , 
I{e)^ f[g-\e)] 



|cosx 



Xel\ ' 



Xel 



(2e— l)cosx sinx — 2-\/e(l — e) |sinx| sinx 



;(e- l)]sin^X- (2e- 1)2 + 4(2e - l)^e(l-e) |sinx|cosx 



etn = g ^(e) = - Ve(l-e)|sinx| + - [1 + (2e - l)cosx] 



1 



q 

a 

2 vr 



Bin G (0, sin ( — 



TTi 



(42) 
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and, 

fg{q) = I fie), q^{0,a) , 
f{e) ^ f[g-He)] ^ 



COSX + Xe2\ 



(2e— l)cosx sinx + 2^/6(1 — e) sin^ x 
y^l + [l + 8e(e- l)]|sinx|siiix- (2e- 1)2 - 4(2e - 1) ^6(1 - e) |sin x| cos x 

ezn = g^\e.) = + Ve(l-e)|sinx| + - [1 + (2e - l)cosx] , 

a 

2 n 
X = ^t, 

e,„ e (sm^f^Vl) . (43) 

The probability P{q < ql) that a perihelion distance q is less than a given value g^ can be 
calculated from 

/•9-L 

^(9 <qL) = / /,(?') dg' , 
Jo 

P{q < ql) = Fg{qL) . (44) 

where Eqs. P^ - ([i5|) are to be used. 



6.1.1. Uniform distribution 

Let us consider, as an example, a uniform distribution of eccentricities. The density and distribution 
functions of perihelion distance are, on the basis of Eqs. (P^. (|^T|) and Eqs. (H^ . 

Jq\Qin) ^ ~ , 

a 

Fqiqrn) = / fqiqL) d<l'in = <iin / a , 

JO 

qtn e (0,a) . (45) 

The distribution function F{qin) fulfills the condition Fq{a) = 1. 
Using Eqs. (|42])-(|43]), we obtain 



/(e) = 


a |cosx + ^ell 

(2e — l)cosx sinx 




X 1 — 


- 2v/e(l-e)|sinx|sinx 




y^l + [1 + 8e(e - 1)] sin^ x - (2e - 






-1)2 + 4(2e- l)^e(l - e) |sinx|cosx 


e = 


a 




X = 


2 TT 

F *' 




&%n G 


(0,sin^(^)), 




and, 






/(e) = 


1 1 

a |C0SX + ^e2| ' 





(46) 



X, 



e2 
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(2e — l)cosx sinx + 2-y/e(l — e) sin^ x 



1 + [l + 8e(e- 1)] |sinx|sinx- (2e- 1)2 - 4(2e - l)\/e(l - e) |sinx| cosx 



a 

2 TT 

X = —t, 



e„, e (sin2 (^1,1). (47) 



As an illustration we take t — P/2. Eqs. (|46| - (|47)) yield fq{q) = 1/a, and, the distribution function 
is Fg{q) = q/a: 

IM) = - , ^ - - 
a 

Fq{q) = - , t = , 
a 

1 P 



fM ^ - , t = 



a 



2 ' 



F,{q) = 1, t=f . (48) 

Eqs. (H^ may be used. 



6.1.2. Distribution function F(e,>,) = e? 



m 



We will consider, as an another example, F(ei„) = e|„ (see Sec. 5.1.2). The corresponding density and 
distribution functions of perihelion distance are, on the basis of Eqs. ([M]) . (PTj) and Eqs. 



/,(*«) = 2 i (l - ^) , 
a V a / 

J^,('Z») = r f^^l'^r.) dqL = *!i (2 - *^ 
Jo a \ a 





%n e (0,a) . (49) 

Using Eqs. (|32])-(|13]), we obtain 

Mq) = -/(e), qe{0,a) , 
a 



^ l + (2e-l)cosx - 2^e(l-e)|sinx| 

^^^^ ICOSX + ^ell 

(2e— l)cosx sinx — 2-\/e(l — e) |sinx| sinx 
y/i + [l + 8e(e- l)]sin2x- (2e-l)2 + 4(2e - 1) v/e(l - e) |sinx| cosx 



e = l - ^ 
a 


: 


27r ^ 




e„ e (0,sin^ 





and, 

/g(g) = - /(e) , qe (0,a) , 
a 



^ l + (2e-l)cosx + 2^e(l-e)|sinx| 
^^^^ Icosx + X,2\ 
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(2e— l)cosx sinx + 2-\/e(l — e) sin^ x 



X. 



e2 



1 + [l + 8e(e- 1)] |sinx|sinx- (2e- 1)2 - 4(2e - 1) ^e(l - e) |sin x| cos x 



e = l-^, 
a 

2 TT 

X = ^t, 
e» e (sin2(^^^,l). (51) 

As an illustration we take t — P/2. Eqs. (l5ni) - (|51l) yield fq{q) = 2 q/a^, and, the distribution 
function is Fq{q) = {q/a)"^: 

Uq) = 2 i (l - ^) , t^Q, 
a \ a/ 



a \ aJ 



a 



a a 2 



PM = (-) ^ t = - 



q e (0, a) . (52) 

The situation represented by Eqs. (|5^ is illustrated in Fig. 7. 



6.2. P{q < qi) for a random time 

Orbital evolution of perihelion distance is given by the evolution of eccentricity, since the relation 
between the semi-major axis a and eccentricity eis<7 = a(l— e), and a is a constant in secular 
evolution due to the tidal effect of the Galaxy. We have, on the basis of Eqs. ([T^ . 

lit) = 2\ + sin (^ IT* - 2 " 
5,n e (- 1, + 1) . (53) 

The probability P{q < q^) that a value q is smaller than a value q^^ in the sense of Sec. 5.2, is 

Pil <qL) = ]: - - arcsin (l - 2 ^) . (54) 

li qL/a <^ 1, then Eq. (|M)) reduces to 

P{q <qL) = - \f^ , qL/a « 1 . (55) 

TT \ a 

If one uses the form g — all — 2|t — tinlimod P)/P\ instead of Eq. (|53l) . then P{q < qi,) — qL/a. 
Eqs. ([51)l - (j55p hold for the case when gravity of the Sun and galactic tide are considered. 

6.3. Distribution of perifieiion distances -galactic tides and results of numerical calculations 

Fig. 8 illustrates distribution of perihelion distances, histogram and cumulative number of comets 
A'^(< q) as a function of a cometary perihelion distance q. The results presented in Fig. 8 correspond 
to the gravity of the Sun and the effect of galactic tide, only. Numerical integrations of equation of 
motion for 4.5 x 10^ years for semi-major axis a™ = 5 x 10^ AU, inclination with respect to the 
galactic equatorial plane iin — 90 degrees and various orientations of the initially almost circular orbits 
(rotation angles - longitudes of the ascending node - 0, 45, 90, ..., 270, 315 degrees) were performed. 
The histogram for the zone g < 30 (40) AU suggests that the (7— distribution corresponds to the 
uniform distribution and it is consistent with the observational data for the group of "new" comets 
presented by Fernandez and Gallardo (1999, Fig. 1). 
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Fig. 7. Density fq{q) and distribution Fq{q) functions of perihelion distance for two instants: t — 
and t = P/2, where P is the oscillation period of eccentricities (perihelion and aphelion distances, 
angular orbital elements). Galactic tide and gravity of the Sun play role. Semi- major axis is normed 
to 1, for simplicity. 



On the basis of the histogram in Fig. 8 one could come to the conclusion that the distribution 
function of perihelia is given by Eq. (1491) and the situation corresponds to the idea of Fernandez 
and Gallardo (1999), Hills (1981) and Jeans (1919), as it was mentioned in Sec. 6.1.2. If the idea is 
physically correct, one should await that the cumulative number of comets iV(< q) is given by the 
formula represented by Eq. P^j) . The least-square method fit to the data yields 

7V(< q) = Aq + B q'^ , q <$: a , 
A = (1.08 ±0.03) AU~^ , 



B = - (4.9 ± 0.3) X 10"^ AU"2 
The comets with a = 5 x 10^ AU are used. Eq. (gH) yields that A = 2 / a and B 



(56) 

' 1 / a^ and 

their ratio is \B\/A = 1 / ( 2 a ) = 1 x 10^^ AU~^ Eq. ([56]) yields \B\/A = 4.54 x 10"^ AU"\ i.e., 
454-times greater than the value corresponding to the idea of Jeans (1919), Hills (1981) and others. 

On the basis of Sec. 6.2 we have tried also the fit of the form iV(< q) = C q''. The least-square 
method yields 

N{<q) 



C-^q"' 



q-^ a 



(3.055 ± 0.335) AU"'' , 
7 = 0.657 ±0.026 . (57) 

The value of the exponent 7 lies between the values 1/2 and 1 discussed in Sec. 6.2. The result suggests 
that 7 equals 2/3. If we use the least-square fit of the form iV(< q) ^ C q^^^ ^ we obtain 

N{< q) = C g2/3 , g < a , 
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C = (2.929 ±0.032) AU"^/^ . 
The error of C is 1%. 



(58) 
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Fig. 8. Histogram and number A''(< q) of comets with perihelia smaller than a value q of the perihelion 
distance. Only cases for the inner part of the Solar System are depicted, q < 100 AU. 

The idea of Fernandez and Gallardo (1999), Hills (1981) and Jeans (1919) (and others) is based 
on the density function /(e) — 2 e or fq{q) — { 2 / a ) ( 1 — q/a ). However, also other functions 
can produce uniform distribution in perihelion distances for small q, e.g., /(e) = ( fc + 1) e*^, fq{q) 
= [il + k) / a]{l-q / a)\ F^iq) = 1- {1- q / a )'=+! (fc > - 1); Fq{q) = (fc + 1 ) g / 
a for small q. Moreover, the distribution treated by Jeans (1919) yields not only a density function 
for eccentricity, but also a density function for semi-major axis. The density function for semi-major 
axis is, according to Jeans (1919), proportional to y/a. This is not consistent with the idea that the 
distribution in semi-major axis is a", where a € (— 4, — 2) (Fernandez and Gallardo 1999). Thus Eq. 
((55|1 may be more physical than Eq. (1551) . 

6.4. Distribution functions Fa, e(a, e), Fa, q{ci,q), Fq{q), Fa{a) 

We want to find distribution functions Fa, e(a, e). Fa, q{a,q) and Fq{q), where a, e and q are the 
semi-major axis, eccentricity and perihelion distance. 



fa{a')da' / f{e')de', 
Jo 



6.4.1. Distribution function Fa, e{o.,e) 

If the density function fa, e{a, e) can be written as the product faici) /(e), then 

Fa, e(a, e) = 

e e (0,1). 
Let {a ^ - I) 

fa{a) - (a + 1) (a^+i - < 
/(e) = (fc + l)e\ fc>0, 

d t \CLmin7 ^riiax) 5 

e G (0,1). 
Eqs. dSni-dSni) yield 



(59) 



(60) 



Fa, e(a, e) 



-,a+l 



a+1 _ a + 1 



a+1 

7nin pk+1 
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e e (0, 1) . (61) 

The case a = — 1 would yield 

/a(«) = (In flmai: " lna„ii„)~ a" -^ , 
i^a, e(a,e) = — e + , 

e e (0, 1) , (62) 

instead of Eqs. (pTj) . 

6.4.2. Distribution function Fa^ q{a,q) 

We have to use e = 1 — q / a in Sec. 6.4.1. 
In general, we have 

Fa, q{a,q) = / / fa, e{ai,e) de dai , 

J a^nin J 1 — q/ a\ 

q e (0,a™„) , 
Fa, q{a,q) ^ / /a, e(ai,e) de dai 

Jamin Jo 

a /.I 



/ / fa, e(ai,e) de dai , 

^g Jl — q/ai 



(X t \G"mini ^max) : 

q e {amin,a) ■ (63) 

If we consider the conventional approach k = 1 (e.g., Hills 1981, Fernandez and Gallardo 1999), 
then Eqs. ^ yield 

Pa, q{a,q) = ^+i" _ „+i i - 9 («" - «"un) - — -T 9^ (a""^ - a™7„^) 

G t \et"rnin 7 ^max / 7 

g e (0,ami„) , 
-Fa, q{a,q) — „+i _ Q+i 






" + 1 f2 . „ „. 1 , . „_i 



'^"^ax ^min ^ ^ a 1 



q e {amin,a) ■ (64) 

If fc is more general, not only the special case fc = 1, then Eqs. (|63p yield (a 7^ — 1) 

^a,,(a,g) = ,+ " + ^+1 r a? j 1 - (1 - -^y^'l d«,= 



fe+i /7, , i\ / 1^^+l 



_ g + l Y^' A + A (- 1)'+^ , . , 



"'min ) I 



g e (0,ami„) , 

O fc \P"rnin7 (^vaax) 7 



20 Klacka et al.: Properties of the Oort cloud 



a+l _ a+1 

y "mm 

a, q\^,H/ - a+1 _ n + 1 
Umax "-min 



Fa., qia,q) 



+ .+r\+i f^^^A'- (i-f)'"^- 

Umax — Umin •' Q { \ "1 / J 

a+1 _ a+1 
y "man 

a+1 _ Q + 1 
Umax "'^^f,,^ 



^max ^min /— i ^ / 

q e {amin,a) , 

fl € \Umin,0'max) ■ (""J 

Of course, fc may not be an integer number. 

The conventional statements on the Oort cloud are: a = — 3/2 or a S ( — 4, — 2 ) and /(e) = 
2 e (e.g., Duncan et al. 1987, Bailey 1983, Fernandez and Ip 1987, Fernandez 1992, Fernandez and 
Gallardo 1999). We have obtained some results on the values of a and Umin already in Sec. 2, see the 
last paragraph of Sec. 2. Also results from Sees. 9 and 10 will be helpful. 

6.4.3. Distribution function Fq{q) 

We have found that Fa. q{a, q) is represented by Eqs. (|M1) if /(e) = 2 e. If we are interested in Fq{q), 
then we have to consider a e {amin,amax)- We obtain Fg{q) = Fa^ q{amax,q)-, or 

1 o / -,,-1 Q-IN 



Fq{q) 


= 


a + 1 f 


2 

a 


Q (a^ax - 


^min 


Q + 1 _ Q + 1 

Umax u-min ^ 


1 


e 


{Q.Umin) , 








Fq{q) 


= 


Q+1 _ Q+1 
y "'min 

Q + 1 _ Q + 1 

Umax '^min 












a + 1 


{ 


2 / Q 


- q" 




1 a+1 Q+1 
«maa; U , 



_ „2 I' Q-1 _ Q-n 



q' («mai - «""'] 



a- 1 

g G {amimUmax) ■ (66) 

Eq. (|66p improves the result of, e.g., Hills (1981), Fernandez and Gallardo (1999). 
The case of arbitrary fc, given by Eqs. (|55|). leads to 

FM) = a+1 ^ ' Q+1 r^^i I 1 - (l - Z-T'] dai = 
Umlix - a„7„ Ja„i„ [ V ai/ J 

^' ^' ^ '^ ^^^^ ?' («m+i-' - a^^-'^ 



a$^ti - a"+i ^ V ^ ; a + 1 - / ^ ^ "'^^ 



imn 



q € (0, a 



Q+1 _ Q+1 
y "miTt 

Q + 1 _ Q + 1 

Umax u^^^^^ 



a + 1 / „ I / q ^ 



1 - 1 - ^ rfai = 



' Q+1 Q+1 , - I \ n 

Umax U^j^^ J a I \ "1 

Q+1 _ Q+1 

y u.^.^.|^ 



a+1 _ Q+1 
Umax U^^^ 



« + l V^A + 1\ (- 1)'^^ J r^a+i-i „"+i-;n 






y V"max y 



q e (amm,amaa;) • (67) 
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As a consequence, Fq{q) ex g^, for small q, holds for arbitrary fc, not only for k ~ 1 or k integer. 
Thus, the result of observations (inner part of the Solar System) yielding Fq{q) ex q^ cannot be used 
as an argument for /(e) = 2 e. 

6.4.4. Distribution function Fa{a) 

The distribution and density functions Fa{a) and /a (a) are 



Fa{a) = 



a+l „"+l 



a"~^" - a 



'TTiin 



Q + 1 _ a+1 ' 
'J'max "'min 



faia) = (a + 1) 



Q + l _ a+1 ' 
Umax "'min 



a G {amin,a,max) , (68) 

if Eqs. (l6Tt and relation Fa{a) — J^ fa{a') da' are used. 

7. Minimal perihelion distances and inclinations for comets with "moderate" initial 
inclinations 

In order to determine some basic properties of secular orbital evolution of a comet under the action 
of the galactic tide we numerically solved the system of Eqs. (13)-(17) in Pastor et al. (2009). We 
consider two values of the semi-major axis, 10 000 AU and 20 000 AU, as examples. We do not consider 
greater values of the semi-major axis since the analytical approach to the secular evolution of orbital 
elements does not yield results equivalent to the results obtained by detailed solution of equation of 
motion (Pastor et al. 2009). 

Table 3 presents numerical solutions qmin and inclinations imin corresponding to the minimal 
perihelion distances qmin- Also oscillation periods P are given in Table 3. By the term "inclination" 
we mean the inclination with respect to the galactic equatorial plane. Table 3 also shows initial 
conditions of the numerical integrations. Initial argument of perihelion is used only from interval Wm 
S (0, 180°). This is is sufficient since orbital evolution for Wj„ = oj + 180° is identical to the evolution 
for ojin = a; if the values of other initial orbital elements are fixed (Pastor et al. 2009). We found that 
evolutions of eccentricity, argument of perihelion and inclination are not very sensitive to initial value 
of the ascending node for the initial inclination ii„ < 80° (under the assumption that the values of 
other initial orbital elements are fixed). Graphical evolutions of the ascending node for various values 
of i^in are approximately parallel, in this case (the parallelism does not hold for the case im « 90°). 
This is the reason why the initial ascending node equals to zero for all numerical integrations in Table 
3. Table 3 also shows that the minimal perihelion distance is a decreasing function of im, if the values 
of other initial orbital elements are fixed. In reality, comets with initial inclinations close to 90° can 
get into the interior region of the Solar System due to the galactic tide. 

Fig. 9 depicts inclination imin at minimal perihelion distance as a function of minimal perihelion 
distance qmin- The values are taken from Table 3. imin is an increasing function of qmin for a given 
ain and iin- Comparison of the values of qmin at Oi„ = 10 000 AU and 20 000 AU, for a fixed Wm and 
«i„, yields 

Qmin 1 ^in 1 /nr\\ 

« . (69) 

Qmin 2 O^in 2 

If ain and i^ are fixed, then the values tUin = w and ujin = 180° — w yield approximately equal values 
of qmin and imin (scc also Table 3). This is the reason why we can see only five distinct points, not 
eight, in Fig. 9. This property is not related only to the values of qmin and imin'- the complete orbital 
evolutions of the comets are very similar for Wj„ — uj and Wj„ — 180° — uj. This property divides the 
interval w™ g (0, 180°) into two parts with different behavior. If we consider only the values uJin G (0, 
90°), then qmin and imin are increasing functions of Wm, if ain and iin are fixed. If ujin £ (90°, 180°), 
then qmin and imin are decreasing functions of u)in - An explanation of this property is probably in the 
antisymmetry of the secular time derivatives after transformations uj—>-Tr — uj,(l—>-Tr — Cl (Pastor 
et al. 2009). Since the orbital evolution is not very sensitive to the initial value of the ascending node, 
we obtain similar orbital evolution for a;.;„ = 180° — uj as for a;.;„ — uj. 
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Table 3. Inclinations imin corresponding to minimal perihelion distances qmin and oscillation periods 
P as a result of numerical integration of the system of differential equations for secular evolution of 
orbital elements. Time of integration is tmt 2 x 10^° years. The initial semi-major axis ai„, eccentricity 
Cim argument of perihelion Wj^, longitude of the ascending node fim and inclination ii„ are shown. 



8. Distribution in the ecliptical inclination 

If we want to consider cometary orbital elements for the inner part of the Solar System, then we have 
to consider initial inclinations with respect to the galactic equatorial plane close to 90 degrees (see 
also Sec. 7). Gravity of the Sun and Galaxy will be considered. 

We used numerical solution of the system of equations for secular evolution of orbital elements 
given by Eqs. (13)-(17) in Pastor et al. (2009). We used the equations to measure minimal perihelion 
distance qmin of a comet at the first three returns of the comet to the inner part of the Solar System. 
At the time when the minimal perihelion distance occurred we measured also the inclination imin 
and the longitude of the ascending node ^min- The results are shown in Table 4. Initial conditions of 
numerical integrations are also shown. Numerical integrations with equal values of initial semi-major 
axes, eccentricities and arguments of perihelion have equal oscillation period. The oscillation period for 
all numerical integrations in Table 4 is P « 1.3 x 10^ years. Change of the longitude of the ascending 
node fi of the comet with initial inclination ii„ k, 90° during the return of the comet to the inner part 
of the Solar System is always ~ ±7r. The last column shows two "stable" values fulfilling the fact that 
f2 lies between the values at perihelion. If the first value is smaller/greater than the second one, then 
ri increases/decreases from the first value to the second one. 
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Fig. 9. Inclination to galactic equatorial plane at minimal perihelion distance Qmin- Square is used for 
a„ = 10 000 AU and z„, = 30°, triangle for ai„ = 10 000 AU and «„ = 60°, star for a™ = 20 000 
AU and ?„ = 30° and circle for «;„ = 20 000 AU and «„ — 60°. The greater «i„, the smaller qmin 
and the greater imim for a given a^- 

We are interested in an inclination i of the cometary orbital plane with respect to a reference plane. 
The reference plane is defined by the inclination iq and the longitude of the ascending node flo with 
respect to the galactic equatorial plane (and a given reference direction) . If the galactic inclination of 
the comet at its perihelion is imin (galactic longitude of the ascending node flmin), then the inclination 
i at perihelion position can be calculated from the following equation: 



cosio cos imm + sin Jo sini„j„j cos{Qo - ^min) 



(70) 



Table 4 shows that the values of inclination with respect to the galactic equator (galactic inclination) 
are always close to 25° or 155°. The number of comets with galactic inclination close to 25° is ap- 
proximately equal to the number of comets with the inclination close to 155°. The values of Qmin are 
practically random (see also Table 4) . These facts yield for a large sample of comets 
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{cosi) 
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1 

2' 
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2 ■ 



sin2 



(71) 



The results state that the density and distribution functions correspond to isotropic distribution, the 
average value of cosi equals zero and the average value of i is tt / 2. The results are consistent with 
Eqs. (©. 

The results presented in Eqs. (I7T|) hold for inclinations with respect to any reference plane ig ^ 0, 
also for the ecliptic plane. The distribution of galactic inclinations for comets in the inner part of the 
Solar System is not isotropic, if only gravity of the Sun and Galaxy are considered: i{galactic) « 25° 
or 155°. 

9. Gravity of the Galaxy, Sun and Jupiter 

This section presents results obtained when also gravity of the planet Jupiter is included. 



9.1. Evolution of orbital elements 

Evolution of semi-major axis and eccentricity is depicted in Fig. 10. The action of the Sun and Galaxy 
produces constant semi-major axis of a comet in the Oort cloud, as for secular orbital evolution. The 
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Table 4. Minimal perihelion distances qmin during the first three returns of a comet to the inner 
part of the Solar System. Inclinations imin and longitudes of the ascending nodes ^min corresponding 
to minimal perihelion distances are given. Initial semi- major axis am, eccentricity ei„, argument of 
perihelion ujin , longitude of the ascending node r2i„ and inclination Zi„ are shown. The longitude of 
the ascending node rapidly changes around r2i„ within the interval presented in the last column. 



presence of Jupiter may cause sudden changes in semi-major axis. Moreover, the planet caused that 
the period of oscillations in eccentricity (perihelion and aphelion distances, inclination, ...) decreased 
to almost one half of the period found without the action of the planet. 



i4it iiiiwiiWi>>i J i>4ii A iiJii 
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Fig. 10. Evolution of semi-major axis and eccentricity for a comet under the gravitational influence 
of the Sun, Galaxy and Jupiter. 
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9.2. Distribution in perihelion distance 

As it is presented in the previous subsection (see Fig. 10), the presence of Jupiter may cause more 
frequent returns of a comet to the inner part of the Solar System than it is in the case when the planet 
is ignored. Moreover, our calculations confirmed the importance of the planet also in another type of 
computational experiment. We considered the same initial conditions (orbital elements - semi-major 
axis 5 X 10'* AU, inclination with respect to the galactic equatorial plane 90 degrees, eccentricity close 
to zero, ...) for the two cases, one without Jupiter and one with the action of Jupiter. We gathered 
long-period comets with perihelion distance q less than 100 AU (two groups of comets do not belong 
to the set of comets: i) comets with q < 0.01 AU, and ii) comets ejected from the Solar System due 
to the close approach to Jupiter). While the ignorance of Jupiter yielded that 5.0% of the comets 
exhibited q < 5 AU, the action of Jupiter, Sun and Galaxy yielded that 28.6% of the comets exhibited 
q < 5 AU. The corresponding percentages for g < 10 AU are: 10.0% of the set without Jupiter and 
33.3% of the set with Jupiter. Thus, our model of the Galaxy relevantly changes the conventional 
result that 'planetary perturbations do not significantly alter the perihelion distance' (see, e.g., Dones 
et al. 2004, p. 161). 

We have already mentioned that the above presented results hold for comets which still belong to 
long-period comets. Although nonnegligible part of comets was ejected from the Solar System due to 
the close encounter with Jupiter, the conventional statement that 'only about 5% of the new comets 
are returned to Oort cloud distances of 10** — 10^ AU' (Weissman 1979; Dones et al. 2004, p. 157) 
is not consistent with our numerical calculations. Our results show that most of the new comets are 
returned to the distances of (lO'' - 10^) AU. 

If we use the least-square fit of the form N{< q) = C q^^^ , where q is the perihelion distance and 
A''(< q) is the number of comets with perihelion distances smaller than the value q, we obtain 

N{< q) ^ C q^'^ , g < 100 AU , 



C = (2.248 ±0.103) AU"^/^ . 
The error of C is 4.6%, greater than the error obtained without the action of Jupiter. 



(72) 
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Fig. 11. Histogram and number A^(< q) of comets with perihelion distances smaller than the value q. 
Only cases for the inner part of the Solar System are depicted, q < 100 AU. 



9.3. Distribution in ecliptical inclination 

Fig. 12 depicts distribution in inclination with respect to the ecliptic. The action of Jupiter is marginal, 
although some difference between Figs. 6 and 8 exist. The isotropic distribution function (see Eqs. 
[S]) of ecliptical inclination - inclination to the ecliptic - is in a good coincidence with the calculated 
data. However, 70% of the calculated orbits exhibit prograde orbits (ecliptical inclination is less than 
90 degrees). 
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Fig. 12. Distribution function of cometary ecliptical inclination when comets are situated at their 
perihelia. Only cases for the inner part of the Solar System are depicted, q < 100 AU. The dotted line 
corresponds to the isotropic distribution. 



One would await that gravity of a planet can influence the ecliptical inclination mainly in the cases 
when smaller perihelion distances occur. If the effect of the planet would not exist, then we should 
await no correlation between the perihelion distance and the ecliptical inclination. Our calculations 
show that the inclusion of Jupiter leads to the coefficient of correlation 



r{q,ieci) = 0.366 ± 0.135 



(73) 



The percentage probability that the correlation coefficient for the used set of data is greater than the 
given value is less than 2%. 

10. Mass of the Oort cloud 

Current estimates of the mass of the Oort cloud of comets are about (3.3 — 7.0) masses of the Earth 
{Me), although the value of 38 Me has also appeared (Dones et al. 2004, p. 162). 



10.1. Observational data and the distribution function Fa{a) 

Using observational data on the original semi- major axes of the long-period comets taken from Marsden 
and William's (1997) catalogue, we have found that the theoretical fit for the density function faia) 
in the exponential form a" holds only for semi-major axes smaller than about 4 x lO'' AU (see Fig. 
13). Theoretical fit for the data, depicted in Fig. 13, yields for the distribution function 



(74) 
0. Eqs. 



Fa{a) = — ) , a < ao , 

aobs + 1 = 1.133± 0.044 AU , 
ao = 4 X lO'' AU . 

Relative error of the exponent is 3.9%. Comparison with Eqs. (|68|) yields a„iax — o-o, a-r, 
(|74[) give that the observed density function of semi-major axis is 



faia) 



^obs 



ao 



aobs 



a-o 



a < ao , 



4 

30 ' 

4 X lO'' AU 



(75) 
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Fig. 13. Dependence of the cumulative number of long-period comets on semi-major axis a [AU] in 
logarithmic scale. The dependence is linear for a < 4 x 10"* AU. 

If we would like to fit the curve in Fig. 13 up to a = 1 x 10^ AU using the approximation given 
by Eqs. dlS]), then 



1 = 0.449 ± 0.049 , 



(76) 



We could fit the curve in Fig. 13 up to a = 1 x 10^ AU using the approximation given by Eqs. 
with greater values of amin- However, large errors of the exponent a -t- 1 exist in these cases. For 
example, the value amin = 0.1 x 10^ AU yields relative error of a -1- 1: about 160%. 

Let us calculate mean value of the semi-major axis. On the basis of Eq. (|68p we can write 



(a) 



a' fa{a') da' 



a 



1 1 



= On 



xP'min / ^max ) 



a+2 



1 y^min/ ^7nax) 



a+1 



(77) 



Figure 14 depicts the value of (a) / a^ax as a function of the exponent a for the cases found by 
various authors: a = — 3/2 (Duncan et al. 1987), a€i( — 4, — 2) (Fernandez and Ip 1987, Fernandez 
1992, Fernandez and Gallardo 1999). Also the value a = — 0.55 given by Eqs. (|75|) is considered. The 
value of (a) / amax for the cases a = — 3/2 and a G ( — 4, — 2 ) are calculated under the assumption 
that amin/amax = 0.1 (see the solid curve and the square in Fig. 14). The case a — — 0.55 considers 
two possibilities, amin/amax = 0.1 (triangle in Fig. 14) and amin ~^ (star in Fig. 14). Our result 
represented by Eqs. ([7S)l yields the value of (a) / amax consistent with the value obtained from the 
model by Duncan et al. (1987). The case a = -I- 1/2 (Jeans 1919) yields more than two times greater 
value than the case a = — 3/2. 

10.2. Number of comets in the Oort cloud 

We will present two different accesses to the estimation of the number of comets in the Oort cloud. The 
first one is based on the comparison between the frequencies of cometary returns followed from the new 
physical model and the standard model of galactic tides. The second access considers the distribution 
of comets in the semi-major axis a when the density function is taken in the form proportional to a". 



10.2.1. The first access 

We have already mentioned the relevance of the physical model of the galactic tide (see Sec. 9.2 and 
the paper by Komar et al. 2009). As a summary, we have that the number of oscillations of orbital 
elements is 10/6 times higher than for the standard model of the galactic tide (see Figs. 2 and 7 in 
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Fig. 14. The ratio of the mean value of semi-major axis to amax as a function of the parameter a 
for various models. If Umin / o-max = 0.1, then the solid line corresponds to the model discussed by 
Fernandez and Gallardo (1999), the square to Duncan et al. (1987), the triangle to a = — 0.55. The 
star corresponds to Eqs. (1761) a = — 0.55 and amin — ^ 0. 



Komar et al. 2009). On the basis of Sec. 9.2 we know that the number of comets coming to distances 
less than 5 AU is 28.6/5.0 times higher than it is in the standard model and the frequency of the 
returns may be even in 15/10 times greater (see Fig. 11). Thus we obtain that the number of occurence 
of the long-period comets in distances less than 5 AU can be 15-times (10/6 x 15/10 x 28.6/5.0 = 
15) greater than the standard model offers (see, e.g., Dones et al. 2004, pp. 161-162). 



10.2.2. The second access 



On the basis of Eqs. p^ and Eqs. ([77]) we obtain that the ratio between the frequencies of cometary 
returns into the inner part of the Solar System for a — — 3.5, Omm / o,max = 1/10 (Dones et al. 2004, 
Fernandez and Gallagher 1999) and a = - 0.55, a,„i„ -> (see Eqs. [751) is /(0.55) / /(3.5) = ((0)3.5 
/ (0)0.55)'^ = 7. On the basis of Sec. 9.2 we know that the number of comets coming to distances less 
than 5 AU is 28.6/5.0 times higher than it is in the standard model and the frequency of the returns 
may be even in 15/10 times greater (see Fig. 11). Thus we obtain that the number of occurence of the 
long-period comets in distances less than 5 AU can be 60-times (7 x 15/10 x 28.6/5.0 = 60) greater 
than the standard model offers (see, e.g., Dones et al. 2004, pp. 161-162). 



10.2.3. Discussion 

The two results were obtained, partially, in two different ways. The first case has not considered any 
density function in semi-major axis a. However, it was based on numerical calculations presented by 
Komar et al. (2009) for a = 5 x 10^ AU, while more correct access is to use the mean value (a) as 
it is in the second case. In any case, both accesses yield values which are more than 10-times smaller 
than the conventional values (see, e.g., Dones et al. 2004). 



10.3. Our estimate of the mass of the Oort cloud 

On the basis of the results discussed above we can come to the conclusion that the mass of the Oort 
cloud is less than 0.5 mass of the Earth (maybe, even 0.1 Me). However, inclusion of perturbation by 
close stars and interstellar clouds may change the result. 

We can present more exact calculation based on the above presented values. The current estimate 
of the mass of the Oort cloud of comets is (3.3 — 7.0) masses of the Earth (Me), (Dones et al. 2004, 
p. 162), or M' = ( 5.1 ± 1.9 ) Me. Our two, partially independent, results yield that the real mass 
M is (15 — 60)— times smaller than the value M', or, M is d = (37 ± 23)— times smaller than M'. 
Thus, M = M' I d. Using also error analysis, we finally obtain M = ( 0.14 ± 0.10 ) Me. 



Klacka et al.: Properties of the Oort cloud 29 

We have not considered the effect of nearby stars (interstehar clouds) . Let the stars can generate 
A^— times higher number of observable comets than the conventional/standard model of galactic tides 
offers. The real number of observable comets, due to the action of more realistic galactic tides, is [ 
(15-60) -|-A'']/(1 + A^) greater than the conventional estimate. Thus, the mass of the Oort cloud is 
( 1 + iV ) / [ (15-60) -I- N ] lower than the conventional estimate. The formula yields, as an example: 
M w (1/20) M' w (1/4) Me for AT = 1, M « (1/13) M' w (1/3) Me iov N = 2, M k (1/7) M' 
« (5/7) Me for TV = 5, M « (1/5) M' w 1 Me for N = 10. If we take into account the result of 
Rickman et al. (2008, Fig. 2), then we should use, as an approximation, iV = 1 (we remind that the 
authors use a = —1.5 which significantly differs from 0). 

11. Summary 

The paper presents various results on the Oort cloud of comets if gravity of the Sun, Galaxy (and 
Jupiter in Sec. 9) are considered. 

Sec. 2 discusses the (density) function of inclination. Besides the first theoretical part, the dominant 
part deals with the inclination with respect to the galactic equatorial plane. The important result shows 
that the exponent a, characterizing distribution of comets in the Oort cloud as a function of semi- 
major axis, should be equal to — 1 if amm corresponds to (10 — 20) x 10'^ AU. This is not consistent 
with the cases treated in the literature (Duncan et al. 1987, Fernandez and Ip 1987, Fernandez 1992, 
Fernandez and Gallardo 1999). If Omm ^ (10 — 20) x 10'^ AU, then various values of a are admitted 
in our treatment of the distribution in inclinations. 

Sees. 3 and 4 present simple relations for some dependencies. They both improve the published 
results and found new results. The relation between the semi- major axis a and oscillation period P is 
some kind of analogy to the third Kepler's law. The relation reads a'^ P = 1 if a and P are measured 
in natural units. 

Sec. 6 seems to be of theoretical character. However, its results are applied to practical problems. 
One question remains open: Why various authors consider distribution in eccentricity in the form 
treated by Jeans (1919) as the relevant theoretical access but the other part of the Jeans distribution, 
corresponding to the density function for semi-major axis, is ignored? Results of Sec. 6 represented by 
Eqs. ([M|) - (|57|) are improvements of the published results (see e.g., Hills 1981, Fernandez and Gallardo 
1999). The result represented by Eq. (p5)) are used in Sec. 10. 

Sec. 6 presents the effect of galactic tide to some observational quantities. The cumulative number 
of comets with perihelion distances is described by Eqs. ([55)1 and ([7^ . which can be generalized to 
the form A^(< q) = N{< qo) {q/qo) , g, go G ( 0, 100 ) AU. This is equivalent to the distribution 
function 

Fgiq) = i-^) , go = 100AU, ge(0,go). (78) 

This result differs from the conventional result based on the Jeans density function /(e) = 2 e. However, 
if the authors do not agree with the Jeans density function of semi- major axis fa{o.) oc ^/a, then they 
cannot use /(e) = 2 e. Correspondingly, the distribution function cannot be of the form Fq{q) = {q 
/ a) (2 — q/a). as the authors state (see the second of Eqs.[S2]). More correctly, if one wants to use 
/(e) = 2 e, he should be able to present an argument in favor of the choice and the argument must 
be independent of the argument presented by Jeans (1919). Moreover, the results of Sec. 6.4.3 show 
that Fq{q) ex q for the inner part of the Solar System holds for any density function of the form /(e) 
= ( fc -I- 1 ) e*^. A: > 0, not only for k — 1. 

Sec. 9 considers also the gravity of Jupiter. The gravitational action of Jupiter significantly infiu- 
ences the distribution in perihelion distance. The real number of long-period comets may be about 
50-times smaller than the number conventionally considered, as it is discussed in Sec. 10. 

We have already mentioned the problem with the values of the exponent a. The exponent charac- 
terizes the distribution of comets in the Oort cloud as a function of semi-major axis. We have found, 
in Sec. 2 for distribution in inclinations, that only the condition amin <^ (10 — 20) x 10^ AU enables 
a ^ — 1. Sec. 10 deals with the value of the a in a different way. Sec. 10 comes to the conclusion 
that a = — 0.55 and amin -^ 0, see Eqs. ([75)1 . Sees. 6.3 and 9.2 yield Fq{q) ex q^^^ if also q greater 
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than several astronomical units are considered. If this holds also for q > amin, then Eqs. (|67p yields 
a — — \ / i. This is, approximately, consistent with the result of Eq. ((76|) . a = — 0.55. Thus, we can 
conclude that a is about — 1/2, more correctly, a £ ( —0.6, — 0.3 ). 

Consideration of the results of the previous sections yields that mass of the Oort cloud of comets 
is ( 0.14 ± 0.10 ) masses of the Earth if the action of the nearby stars (galactic clouds) is negligible. 
The mass of the Oort cloud is less than 1 Me even if the effect of the nearby stars is important. 

12. Conclusion 

The main results are: 

1. Theoretical number of long-period comets with perihelion distance q < 5 AU is about 50-times 
greater than the conventional approach yields. (Gravity of Jupiter was taken into account in finding 
this result.) Mass of the Oort cloud of comets is, probably, about 1/4 mass of the Earth. 

2. Semi- major axis a and period of oscillations P of eccentricity (and other orbital elements) are 
related as a"^ P = 1 in natural units for a moving Solar System in the Galaxy. The natural unit for 
time is the orbital period of the Solar System revolution around the galactic center and the natural 
unit for measuring the semi-major axis is its maximum value for the half-radius of the Solar System 
corresponding to the half-radius of the Oort cloud. The relation holds for the cases when comets 
approach the inner part of the Solar System, e.g., perihelion distances are less than w 100 AU. 

3. The minimum value of semi-major axis for the Oort cloud is a„ji„ <C 1 x 10"* AU. This condition was 
obtained both from the numerical results on cometary evolution under the action of the galactic tides 
and from the observational distribution of long-period comets. If the density function of semi-major 
axis is approximated by proportionality a", then a is, approximately, — 1/2. 

4. The magnitude of the change in perihelion distance per orbit, Ag, of a comet due to galactic tides 
is a strong function of semi-major axis a, proportional to a®'^^. 

The usage of the density function of eccentricity in the form /(e) = 2 e cannot be argumented 
to be the result obtained by Jeans (1919). Consequence of the Jeans calculations is marginal density 
function of semi-major axis proportional to ^/a and this is not consistent with the conventionally used 
distributions. Moreover, Figs. 7 and 8 in Komar et al. (2009) show that gravitational tides alone would 
produce only values e < 0.9 for galactic inclinations i less than 90 degrees and practically e < 0.8 
for i w 90 degrees: only for very short time intervals e is greater than 0.8 (e is close to 1 during the 
cometary visit of the inner part of the Solar System). 

Using the observations yielding the distribution function of perihelion distance Fq(q) ex q (inner 
part of the Solar System), we cannot come to the conclusion that /(e) = 2 e. Any density function of 
the form f{e) — { k + I ) e'' , k > is also consistent with the observational result Fq{q) ex q. 

Observational data on distribution in semi-major axis a, theoretical analytical results on the dis- 
tribution functions Fa{a) and of a perihelion distance, and, detailed numerical calculations on gravi- 
tational action of the Galaxy lead to a conclusion on Fa{a). The conclusion is that approximation of 
Fa{a) by the form proportional to a"+^ is consistent with all the mentioned methods if a is about — 

1/2. 

The obtained results may be improved to be more accurate. This will require also more robust 
and more detailed numerical calculations (and, also, several planets, not only Jupiter, have to be 
considered). 

Appendix A: Inclination with respect to the ecliptical plane 

(Reference to equation of number (j) of this appendix is denoted as Eq. (A j). Reference to equation 
of number (i) of the main text is denoted as Eq. (i).) 

We are interested in the inclination with respect to the ecliptical plane if we know inclination (and 
longitude of the ascending node) with respect to the galactic equatorial plane. 

The transformation will be found in several steps. 

We will use Xi, i = 1, 2, 3 (and also primed quantities for other reference frames, all with origin 
in the Sun) for coordinate right-handed axes. Summation convention is adopted, i.e. summation over 
repeated indices is assumed. 
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At first, wc will find unit vector normal to the (osculating) orbital plane of the comet. Let the 
orbital plane is characterized with the ascending node fl and galactic inclination i in the unprimed 
system. In order to obtain coordinates in the system where preferred plane is given by the orbital 
plane of the comet, we will make two transformations: 

/ 4 

All = + cosf] , Ai2 = + sinf] , A13 = , 
A21 = - sin fi , A22 = + cosfl , A23 = , 

A31 = , A32 = , A33 = +l , (1) 

and, 

x" - B x' 

Bii = + 1 , Bi2 = Q , Si3 = , 

-B21 = , B22 = + cosi , B23 = + sini , 

B31 = , B32 = - sini , _B33 = + cosi . (2) 

Eqs. (Al) and (A2) yield {dj = Bik At,) 

Cii = + cosfJ , Ci2 = + sinfi , Ci3 = , 

C21 == — sinO cosi , C22 = + cosil cosi , C23 = + sini , 

C31 = + sinf7 sini , C32 = — cos f7 sini , C33 = + cosi . (3) 

The inverse transformation is (transformation matrix C'^ is inverse, in our case transpose, to the 
matrix C): 



x^ 


= 


l^ij Xj , 








= 


+ cosfJ , 


C12 ~ — sinfl cosi , Cj^ = 


= + sin f2 sin i , 


*-^21 


~ 


+ sinO , 


C22 = + cos fl cos i , C23 = 


= — cosO sini , 


'^31 


= 


0, C3^2 = 


= + sin i , C33 = + cos i . 





(4) 
The unit vector normal to the orbital plane of the comet is characterized by the condition 

(x'/,4',4r - (0,0,1)^. (5) 

Its galactic coordinates are, using Eqs. (A4), 

En = (a;i,a:;2,a;3)^ , 

Xi = + sinf7 sini , 

X2 — ~ cos fi sin i , 

X3 — + cosi . (6) 

The top of the vector e^r is a point. The point has the following galactic coordinates 

En = (a;i,a:;2,a;3)^ , 

xi — + cos I cos b , 

X2 ^ + sini cos 5 , 

0:3 = + sin5 , (7) 

where I is the galactic longitude and b is the galactic latitude. 
Comparison of Eqs. (A6)-(A7) yields 

cos I cos b = + sin Q sin i , 
sin I cos 6 = — cos il sin i , 

sin& = + cosi . (8) 
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If n and i arc given, then Eqs. (A8) offer 

b = arcsin(cos i) , 

cos / = + sin fl sin i , 

coso 

sin / = — cos Q sin i . 

coso 

(9) 
We can easily find I from Eqs. (A9). It is sufficient to use the following prescription: 

cos * = C , 

sin * = S , 

S >0 ^ * = arccos C , 

S <0 ^ * = 27r - arccos C . (10) 

We can summarize: We have i and O of a cometary orbit with respect to the galactic equatorial 
plane. The corresponding unit vector normal to the cometary orbit e^r is given by Eqs. (A6). Galactic 
coordinates I and b of the unit vector are given by Eqs. (A7) and they can be found using Eqs. 
(A8)-(A10). 

Equatorial coordinates of the vector e^r are characterized by the right ascension a and the declina- 
tion S: (en eq.c.)i — cosa cos (5, (e^r eq.c.)2 = sina cos (5, (e^v eq.c.)3 = sin (5. Transformations between 
galactic and equatorial coordinates yield 

sin (5 = sin 6 sin 6o + cos 5 cos Sq sin(^ — ^q) j 

cos(q; — Qfo) cos (5 = sin 5 cos Sq — cos 6 sin Sq sin(^ — ^q) ) 

sin(a — Qfo) cos 6 — cos 6 cos{l — Iq) , 

ao == 192.86° , 

6o = 27.13° , 

lo = 33.93° . (11) 

If n and i are given, then Eqs. (A8)-(A11) enable to find a and S. li i ^ 0, then S = Sq and a = 

ao- 

Ecliptical coordinates of the vector e^ are characterized by the ecliptical longitude A and 
the ecliptical latitude /3: {ejs[ eci.c.)i = cosA cos/3, (e^r ec/.c.)2 = sinA cos/3, (e^v ec;.c.)3 = sin/3. 
Transformations between equatorial and ecliptical coordinates yield 

sin /3 — sin S cos e — cos 6 sin e sin a , 

cos A cos/3 = cos (5 cosa , 

sin A cos/3 = sin (5 sine + cos 6 cose sina , 

e = 23.5° . (12) 

If n and i are given, then Eqs. (A8)-(A11) enable to find a and 5. li i ^ 0, then S = Sq and a = 
ao. If a and 5 are known, then Eqs. (AlO) and (A12) enable to find A and /3. If (5 = 90°, then A = 90° 
and /3 = arcsin(cose). 

Finally, the longitude of the ascending node fled and the inclination ied of the cometary orbit 
measured in the ecliptical coordinate system are given as follows: 

+ sin r^ec; s'mieci — cos A cos/?, 
— cos Qeci sin ieci ~ sin A cos /3 , 

+ cos Zee/ — sin/3 . (13) 

The solution of Eqs. (A13) is 

ieci — arccos(sin/3) , 
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siiii7eci = H COS A cos/3, 

Sm led 

cosi7ec( — — sin A cos/3 , (14) 

sm led 

where also Eqs. (AlO) must be used. 

Appendix B: Motion of Jupiter in the ecliptical plane 

(Reference to equation of number (j) of this appendix is denoted as Eq. (B j). Reference to equation 
of number (i) of the main text is denoted as Eq. (i).) 

If we want to take into account also gravity of a planet, we need to find its position in the galactic 
coordinates. Let us consider Jupiter moving in a circular orbit in the ecliptical plane. Using results of 
the Appendix A, we can write for Jupiter coordinates 

x'i J "^ o-J cos(cj,/ t + (^o) , 
X2 J ^ aj sm{ujj t + Lpo) , 

2 TT _i 

a/ 
aj = 5.203 AU , (1) 

where t is the time and ipo is the arbitrary initial phase. The unit vector normal to the orbital plane 
of the comet is characterized by the condition 

i^'lj,4j,4jf - (0,0, If. (2) 

As for the galactic coordinates, we have (see Eqs. A6) 

ieN)gal = (a;i,X2,X3)^ , 

xi = + sin 17,7 sinij , 
X2 = — cosflj sinij , 
X3 = + cosij . (3) 

Eqs. (A4) yield 



Xl J 


= Xl J cosfij — X2 J sinfij coaij 


X2 J 


— x'l J sinfij + x'2 J cosfij cos? J 


X3 J 


= X2 J sinij . 



(4) 

We need the values of Q.j and ij. 

We have {eM)eci = ( 0, 0, 1 )'^ . Transformations inverse to Eqs. (A12) 

sin (5 = sin/3 cos e + cos/3 sin e sin A , 

cos a cos (5 = cos/3 cos A , 

sin a cos (5 = — sin/3 sine + cos/3 cose sin A , 

e ^ 23.5° , (5) 

yield for /3 = 90° 

sin(5j = + cos e , 

cos a J cos (5j = , 

sin a J cos 5 J = — sine , 

c = 23.5° . (6) 
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Solution of Eqs. (B6) is 

dj = arcsin(cos e) , 
aj = 270° , 
e = 23.5° . (7) 

Thus, we have 

(eAr)eg — {cos a J cos(5j,sinaj cos(5j, sin(5j) 
— (0,— sin£, + cose)^ , 
e = 23.5° . (8) 

As for transformation to the galactic coordinate system, we have to use transformations inverse to 
those represented by Eqs. (All): 

sinfe = sin (5 sin Sq + cos 5 cos Sq cos(q! — ao) , 

cos(/ — Zo) cos 6 = cos (5 sin(Q: — ao) ; 

sin(/ — lo) cos 6 — sinS cos Sq — cos 6 sin Sq cos(a — ao) , 

ao == 192.86° , 

So = 27.13° , 

lo = 33.93° , (9) 

which yield, together with Eqs. (B7), 

sin 6,7 = sin 5q cose — cos 6o sine sinap , 
cos{Ij — Iq) cos bj = — cosao sine , 
sin(^j — Zq) cosbj = + cos Sq cose + sin Sq sine sinao . (10) 

Finally, we can write 

{eN)gai = {coslj cosbj,smlj cos 6,7, sin 5j)^ 

. = (sinOj sinij,— cosf2j sinij,cosij) . (11) 

Eqs. (BlO)-(Bll) give 

ij = arccos {sin Sq cose — cos 5q sine sinao} , 

r . . . 2l^/^ 

cos 6.7 = -^ 1 — (sin (5o cose — cos Sq sine sinao) r 

cosilj — Iq) ~ ;— (— cosao sine) , 

cos 0,7 

sin(^j — Zo) = -, — (cos (^0 cose + sin (5o sine sinao) , 

cos 0,7 

sin f2 J ~ H — ; cos Ij cos bj , 

sini,7 

cosfij = : — — sinZj cos6j . (12) 



sinij 



Eqs. (AlO) must be used in finding Ij and ilj. 
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